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Pair production of J/ip as a. probe of double parton scattering at LHCb 
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''Institute for Theoretical Particle Physics and Cosmology, 
RWTH Aachen University D-52056 Aachen, Cermany 

We argue that the recent LHCb observation of J/^-pair production indicates a significant contri- 
bution from double parton scattering, in addition to the standard single parton scattering compo- 
nent. We propose a method to measure the double parton scattering at LHCb using leptonic final 
states from the decay of two prompt J/xp mesons. 
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Introduction. The Large Hadron Collider (LHC) 
provides a unique environment for precise measurements 
of hitherto poorly understood phenomena. Since the flux 
of incoming partons increases with the collision energy, 
there is a high probability at the LHC of multiparton 
scattering, i.e. scattering of more than one pair of par- 
tons in the same hadron-hadron collision. The parton- 
parton correlations and distributions of multiple partons 
within a proton are difficult to address within the frame- 
work of perturbative QCD. Therefore detailed experi- 
mental studies of multi-parton interactions are of great 
importance. In particular, it is widely expected that mea- 
surements of double parton scattering (DPS) processes 
with final states carrying relatively large transverse mo- 
mentum (pt) will provide relevant information on the 
nature of multiple scattering. Probing DPS processes us- 
ing leptonic final states has been discussed in [2] . In this 
Letter, we discuss how observing four-muon final states 
from pair production of J/ip could provide additional ex- 
perimental input. 

Taking advantage of high jet production rate at hadron 
colliders, DPS searches have been performed by the 
AFS UA2 [i], CDF [| and DO Q Collaborations 
in the four-jets (4j) and 7 -I- 3j channels. At the LHC, 
the pair production of muons from single J/ip produc- 
tion benefits from a large cross section. This implies 
a significant DPS production rate of two J/tp particles, 
which can subsequently decay into four muons, leading 
to much cleaner signals compared with those in the jet- 
based studies. In addition, the measurement of these J/ip 
pairs provides complementary information on parton- 
parton correlations, as the hard processes can be initi- 
ated by partons different from those leading to Aj and 
7 4- 3j events. 

A pair of J/ip mesons can also be produced in a sin- 
gle parton scattering (SPS) process, e.g. gg — ?► J/iJjJ/iJj. 
Studies of these SPS processes are expected to provide 
important insights for improving the theoretical descrip- 
tion of single quarkonium production. Understanding the 
DPS contribution is thus an important task in this con- 
text. However, here we are primarily interested in the 
DPS as a signal process and will regard the SPS produc- 
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FIG. 1: Invariant mass distribution of the JJ^ pair 
i'm-j/^j/if,)- The LHCb data are read off the plot in From 
top to bottom, the histograms are the DPS-I-SPS, DPS only 
and three (practically degenerate) SPS predictions including 
different radiation effects. 



tion as an irreducible background. 

Invariant mass distribution of J/V' pairs at 

LHCb. The low pT muon trigger at LHCb provides 
an ideal laboratory for studying four-muon final states 
from the decay of two J/ip mesons. The LHCb Collab- 
oration has recently reported a first measurement Q of 
this process, and compared the two-J/-0 invariant mass 
(wj/i/'J/i/j) distribution with the theoretical prediction 
of [7| from direct SPS production of two J/ ip particles us- 
ing a leading-order (LO) color-singlet result first derived 
in Q . In Fig. [1] we show a comparison between the data 
and an SPS prediction for the process gg — >■ J/'ipJ/ip, 
which we implement in the Monte Carlo event gener- 
ator Herwig++ v2.4.2 0. This process is also imple- 
mented in the event generator DJpsiFDC llCl . We neglect 
the known LO color-octet contributions |llUl2l as they 
are predicted to be negligible at the LHC In the 

calculations, the wave function at the origin takes the 
value |i?(0)p = 0.92 GeV^ We use the MSTW 2008 
NLO PDFs and set the renormalization {^r) and 
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factorization scales (^f) equal to ^m?j^^-\- p^, where 

~ 3.096 GcV is the physical J/?/) mass. The J/ip 
mesons arc produced on-shell and decay isotropically 
into ii'^ pairs with the branching ratio BR{J/il) — !> 
fX^ jjC) = 0.05935. This approximation, ignoring possi- 
ble J polarization effects, is justified if one assumes 
small values of J/ip transverse momentum In our 

studies we require, in correspondence with the experi- 
mental analysis, that the J/ij: mesons have rapidities 
2 < yj/^, < 4.5 and transverse momenta Pt,j/iI> < 10 
GeV. 

The framework of Herwig++ allows us to include ef- 
fects of the initial state radiation (ISR) and intrinsic px 
of initial state partons. We set the root mean square in- 
trinsic pt of a Gaussian model in Herwig++ to 2 GeV. In 
the figures, they are referred to as "shower" and "intrin- 
sic" respectively. As can be seen in Fig. [U the impact 
of these effects on the predicted tuj/^j/^ distribution is 
negligible. 

Fig. [T] shows that the rrij/^iji^ shape of the standard 
SPS contribution does not match the data very well. In 
particular, it peaks at too low a mj/^j/jp value. The 
same conclusion has also been derived in [6[ for the theo- 
retical predictions of 0] . In addition to the SPS contribu- 
tion, in Fig. [1] we also show the invariant mass distribu- 
tions for the DPS process, which will be discussed shortly. 
The DPS contribution has a broader shape and peaks at 
higher values of tjj/^j/^,, leading to much better agree- 
ment between the combined SPS+DPS prediction and 
data than with the SPS contribution alone. The LHCb 
data could therefore indicate the presence of a significant 
DPS contribution to the double J /i}} production process. 

Our predictions for the DPS double J production 
are obtained using an approximation in which the DPS 
cross section (cr'^'p^) factorizcs into a product of two SPS 
cross sections (c^gpg): 



da 



2J/ip 
DPS 



2(Tnff 



(1) 



where 



^^SPS = X! fa.i^o.,f^F)fb{xb, fJ-p) f^^SPS dXadXb , (2) 



is the SPS cross section of a single J/ip-, and 

d^SPS ^ X! OS \^Iab^j/ip+x\'^ dPSj/^+x (3) 



is the corresponding parton level process, with s the 
partonic centre~of-mass energy. The approximation as- 
sumes factorization between the longitudinal and the 
transverse components of the generalized double parton 
distributions and the assumption of no longitudinal mo- 
mentum correlations between the partons in the same 



hadron. In this framework, all the information on the 
transverse structure of the proton is captured in the fac- 
tor (Toff. The assumption of factorization between the 
longitudinal and transverse components of double par- 
ton distributions appears to be not strictly valid within 
QCD [isl - flTj . However, given the small values of longitu- 
dinal momentum fraction x probed in the J/ip production 
at the LHC, this should be a reasonable approximation. 
For the same reason, the approximation of no longitudi- 
nal momentum correlations may also be justified. 

The state-of-the-art theoretical description of direct 
single J/^p production involves the non-relativistic QCD 
factorization approach [18| . Because of the LHCb ability 
to trigger on low px muons, down to 1 GeV, it is essential 
for our analysis to be able to describe the production of 
J/tp with low Pt accurately. Since the fixed-order cal- 
culations fail to provide such a description in this 
regime, we resort to modeling the pt distribution of the 
single J/ip- To do this, we approximate the matrix ele- 
ment for the inclusive production of a prompt J/ip, as- 
sumed to be dominated by the gluon-gluon channel, with 
a crystal ball function of the form 



\Mgg^J/^, + x\'^ — 
2 



PT < (pt) 

ifexp(-^^)(l + ^P^^^)-^ pr > (pt) , 



(4) 



■r/i' 



r/i, 



where K = X Ks/rrij^^. The values of the coefficients 
K, A, n and (pt) are obtained throu gh a combined fit of 

and 



ATLAS 



211, CMS 



da/dpT to the LHCb | 

CDF m data, using MSTW 2008 NLO PDFs with the 
factorization scale fip = + Pt ^"^^ ^^'^ physical 

mass mj/^ = 3.096 GeV. The resulting fit gives k, ~ 
0.6 and A 0.327 for n = 2 and {pr) = 4.5 GeV. In 
particular, ~ 28 for the fit to 55 LHCb data points. 

The parametrization (jl]) is then used to obtain the 
cross section for the production of two J/tp mesons 
through DPS. The corresponding matrix element is im- 
plemented into the framework of Herwig++ , with parton 
shower and intrinsic p^-broadening switched off as the 
fitted data already account for these effects. 

The value of cTcff in the factorized approach to DPS, c/. 
Eq. ([T|), which could be energy- and process- dependent, 
is one of the properties of the DPS that requires a more 
precise experimental measurement. In the calculations 
we use CToff = 14.5 mb, a value obtained by the CDF 
7 4- 3j study . As clearly seen in Fig. [1] double J/ ip 
production at LHCb offers a promising opportunity to 
probe the DPS component. In the following, we propose 
a method to separate the DPS and SPS contributions. 

Measurement of DPS at LHCb. The majority of 
kinematic variables used in the literature to distinguish 
a DPS signal from the SPS background, where all four 
final states originate from a single parton-parton hard 
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FIG. 2: DPS and SPS cross sections for pp J/ipJ/ip + 
X, as a function of minimum px of J/V', constructed from 
(upper panel), and the fraction of DPS events using 
the "shower + intrinsic" SPS calculation (lower panel). Basic 
selection cuts applied. 



scattering, are based on the idea of pair-wise balancing. 
For DPS, the two final state particles from the same hard 
scatter will balance against each other on the plane trans- 
verse to the collision axis, resulting in equal but opposite 
transverse momenta. The balancing is exact at leading 
order, but, as we will see in the following, radiation effects 
can have significant impact and reduce the effectiveness 
of these variables. 

In our simulation, we consider events with four muons, 
each required to have pr > I GeV, and to lie within pseu- 
dorapidity range 2 < 77 < 4.5. Out of the two combina- 
toric ways to form two J/tp candidates, the combination 
with invariant mass closest to the physical J/ip mass is 
chosen. In the rest of this Letter we refer to this set of 
cuts as the 'basic cuts'. Once the muon candidates fulfill 
these cuts, 100% reconstruction and detection efficiency 
is assumed. The same parameter values and PDFs as in 
the calculation of Fig. [1] are used. 

Fig.[2]shows the DPS and SPS cross sections as a func- 
tion of the minimum px of a muon pair, after applying 
the basic cuts. In the upper panel, the effects of the par- 
ton shower and the intrinsic py-broadening on the SPS 
predictions is shown. These effects are significant, due 
to the low invariant mass of the two-J /ip system. In the 
lower panel, we see that the DPS fraction increases as 
minimum pT decreases. Conversely, a cleaner SPS sam- 
ple might be obtained by imposing a cut on the minimum 
Pt of the J/ip- 

The impact of the ISR and the intrinsic py-broadening 
is also demonstrated in Fig.[3l where we show the distri- 
bution of A(/) = A0(/.(+/i^, /i+/i^), the azimuthal angular 
separation between the two reconstructed J/tp^s. As ex- 
pected, the signal distribution is flat, a reflection of the 



FIG. 3: Angular separation A<j} on the plane transverse to the 
beam axis between the two J/ip's, reconstructed from fi'^ ii~ 
pairs, for DPS and SPS contributions to pp — > J/ipJ/ip + X. 
Basic selection cuts applied. 



BR^ X cross sections [pb] at 7 TeV LHCb 




DPS 


SPS 


basic cuts 


2.18 


1.27 


basic cuts -|-|Aj/[ > 1.0 


0.48 


0.07 


basic cuts +|Aj/| > 1.5 


0.11 


0.005 


basic cuts + Ay > 2.0 


0.007 


< 0.001 



TABLE I: DPS and SPS cross sections, including branching 
ratios into muons. The first row shows the cross section af- 
ter basic selection cuts, while the rest shows the impact of 
an additional |Aj/|min cut. For clarity, only the SPS result 
with parton shower and intrinsic pr-broadening turned on is 
shown. 

independent scattering hypothesis. For the background, 
while A</) = TT at the parton level, the distribution is heav- 
ily distorted in the presence of ISR and pT^broadening. 
In particular, the ISR leads to distributions that are flat 
or even peaked towards Acf> = 0. We conclude that vari- 
ables based on pair-wise balancing might not be the best 
tools to distinguish between DPS and SPS in this partic- 
ular analysis. 

However, as wc now demonstrate, it is possible to use 
correlations along the longitudinal direction between the 
two J/ip mesons to extract the DPS signal. The idea re- 
lies on the observation that in order to minimize the in- 
variant mass of the J/tp pair, the SPS background should 
on average be characterized by a small rapidity separa- 
tion (Ay). To see this, note that in a frame where the 
Pt of the JZ-^-pair is zero, 

mj/^j/^ = 1^rr^~^~V^ cosh ^ (5) 

hence a small |Ay| is preferred. However, this constraint 
does not apply to the DPS signal, which implies a broader 
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distribution. The small invariant mass of the system en- 
sures that overall momentum conservation has negligible 
impact on the y, and hence the |Ay|, distributions. 
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FIG. 4: Rapidity separation |Ay| between the two recon- 
structed J/il^'s, with basic selection cuts applied. 
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FIG. 5: DPS and SPS cross sections for J/tp-pair production 
with basic selection cuts applied, as a function of a lower cut 
|Ay|inin on rapidity separation. 

As shown in Fig. |4] the difference in |Ay| distributions 
persists in the laboratory frame. As expected, the DPS 
signal is broader and extends to higher values of \Ay\. 
The distributions are also more stable against radiation 
and intrinsic px effects when compared with the Acjj dis- 
tributions, making the predictions more robust. 

To extract the DPS signal, we apply a lower cut on 
the rapidity separation, \Ay\ > |A?/|min. The variation 
of the cross section with jAj/lmin is displayed in Fig.[5l 
Clearly, the event sample becomes more dominated by 
DPS contributions for higher values of |Aj/|min. A sum- 
mary of the results is displayed in TableUl In the current 



(7 TeV) LHC run, an integrated luminosity of a few fb~^ 
is expected. By selecting the four-muon signal sample 
using the basic cuts, a signal to background ratio of a 
few to one may be achieved. Hence we conclude that 
DPS can be measured at the LHCb in the four-muon 
events already at this stage of LHC running. 

Summary. Precise measurement of double parton 
scattering processes at the LHC is an important step to- 
wards understanding multiple interactions in hadron col- 
lisions. The characteristics of the LHCb detector make it 
particularly well suited to study DPS in the production 
of a J/'0-pair decaying into four muons. We observe that 
the first LHCb data on the invariant mass distribution of 
the J/ '0-pair system might already indicate a significant 
contribution from the DPS production mechanism. The 
studies presented in this Letter show that it is possible 
to measure the DPS component in the four-muon events 
at the LHCb already in the early stages of the LHC run- 
ning, in particular with the help of the proposed rapidity 
separation variable Ay. 
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